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ABSTRACT 

We look for observational signatures that could discriminate between Newtonian and modified Newtonian (MOND) dynamics in the 
Milky Way, in view of the advent of large astrometric and spectroscopic surveys. Indeed, a typical signature of MOND is an apparent 
disk of "phantom" dark matter, which is uniquely correlated with the visible disk-density distribution. Due to this phantom dark 
disk, Newtonian models with a spherical halo have different signatures from MOND models close to the Galactic plane. The models 
can thus be differentiated by measuring dynamically (within Newtonian dynamics) the disk surface density at the solar radius, the 
radial mass gradient within the disk, or the velocity ellipsoid tilt angle above the Galactic plane. Using the most realistic possible 
baryonic mass model for the Milky Way, we predict that, if MOND applies, the local surface density measured by a Newtonist will 
be approximately 78 M Q /pc 2 within 1.1 kpc of the Galactic plane, the dynamically measured disk scale-length will be enhanced by a 
factor of 1.25 with respect to the visible disk scale-length, and the local vertical tilt of the velocity ellipsoid at 1 kpc above the plane 
will be approximately 6 degrees. None of these tests can be conclusive for the present-day accuracy of Milky Way data, but they will 
be of prime interest with the advent of large surveys such as GAIA. 
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1. Introduction 

Our understanding of Galactic stellar populations and kinemat- 
ics achieves regular progress with the introduction of new large 
surveys of stars with photometry, distances, radial velocities or 
proper motions, enabling us to probe the 6-dimensional space of 
stellar positions and velocities (e.g., Perryman et al. 1997; Hog 
et al. 120001 Nordstrom et al. 120041; Famaey et al. 120051 Zwitter 
etal.|2008]l. 

For instance, the shape of the Galactic 3-dimensional po- 
tential has been probed by, e.g., the orbits of the Sagittarius 
stream (Ibata etal 120011 Helmi 120041 Read & Moore 120051 
Haghi et al. 120061 1. the tidal tails of Palom ar 5 (Grillmair & 
Dionatos 2006; Grillmair & Johnson 2006), or the kinematics 
of halo stars. In the solar neighbourhood, the potential can also 
be analysed by measuring the Galactic escape speed of high ve- 
locity stars (e.g., Smith et al. 2007 ; Famaey, Bruneton & Zhao 
120071 ). the force perpendicular to the Galactic plane (e.g., Oort 
[19601 Creze et al Q M1 Siebert et al. 120031 Nipoti et al. 120071 
Holmberg & Flynn 2000 , 2004 ), the coupling between the three 
components of the velocity distribution in the solar neighbour- 
hood (e.g., Bienayme 1999), or the orientation of the velocity 
ellipsoid above the Galactic plane (Qllongren |19621 Hori & Liu 
[19631 Lynden-Bell [l962l Siebert et al. 121)081 . In the future, much 
progress is still expected to be made in our understanding of the 
Galactic potential with the advent of the JASMINE and GAIA 
missions. 

This mapping of the Galactic potential is of great impor- 
tance because it holds the key to our understanding of the be- 
haviour of dark matter on galactic scales. Nowadays, the dom- 
inant paradigm is that dark matter is made of non-baryonic 
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weakly-interacting massive particles, called "cold dark matter" 
(CDM). However, on galactic scales, the observations appear to 
be at variance with a sizeable list of CDM predictions (see, e.g., 
Famaey et al. 2007 ). The observed conspiracy between the mass 
profiles of baryonic matter and dark matter at all radii in spiral 
galaxies (e.g., Famaey et al. 2007 ) rather lends support to modi- 
fied Newtonian dynamics (MOND, Milgrom 1983 ), a paradigm 
postulating that for accelerations below oq « 10~ 10 ms~ 2 the true 
gravitational attraction approaches (g^cio) 1 / 2 , where is the 
usual Newtonian gravitational field. Without resorting to CDM, 
this paradigm is known to reproduce galaxy scaling relations, as 
well as the rotation curves of individual galaxies ranging over 
five orders of magnitude in mass (e.g., Sanders & McGaugh 
2002 ). In particular, the kinematic analysis of tidal dwarf galax- 
ies by Bournaud et al. (2007) is difficult to explain within the 
classical CDM framework, while it is in accordance with MOND 
(Gentile et al. 2007). On the other hand, some problems arise 
with this paradigm in the subgalactic and extragalactic scales 
(e.g. Zhao 120531 Angus et al.l2TjQ71 see e.g., Zhao 120081 Angus 
120091 and Bruneton et al. 2009 for possible solutions). 

The present study examines possible observational signa- 
tures of MOND gravity in the Milky Way, which could be de- 
tected (or rejected) with the advent of large Galactic surveys. 
This is in direct continuity with the works of Famaey & Binney 
d2005l l; Nipoti et al. d200Tb ; Wu et al. J20081 , and McGaugh 
(120081 1 . We build th ese pr edictions using the MOND Milky Way 
model of Wu et al. (2008), and refer the reader to this paper for a 
full description of the modefl This model is based on one of the 



1 W e use the model labelled "MOND g cxt = 0.1a " in Wu et al. 
( 2008 ), meaning that the modulus of the external gravitational field act- 
ing on the Milky Way is chosen to be oqI 100. 
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most realistic possible baryonic mass models of the Milky Way, 
the Besancon model (Robin et al. 12003b . 

Once the gravitational potential of the model is known, one 
can apply the Newtonian Poisson equation to recover the den- 
sity distribution that would have yielded this potential within 
Newtonian dynamics. In this context, MOND predicts a disk of 
"phantom" dark matter (Sect. 2), allowing us to differentiate the 
MOND prediction from these of a Newtonian model with a dark 
halo by (i) measuring dynamically (within Newtonian dynam- 
ics) the disk surface density at the solar radius (Sect. 3), (ii) 
measuring the radial dynamical mass gradient within the disk 
(Sect. 4), or (iii) measuring the velocity-ellipsoid tilt angle above 
the Galactic plane (Sect. 5). We show that none of these tests 
are yet conclusive with present-day accuracy, but they are ex- 
tremely useful to future high-precision astrometric and spectro- 
scopic surveys. 

2. The "phantom" disk of MOND 

Bekenstein & Milgrom (1984) proposed a modification of the 
Newtonian dynamics to reproduce the flat rotation curve of disk 
galaxies without the introduction of enormous amounts of dark 
matter on galactic scales. They modified the Poisson equation 
by relating the gravitational potential cD to the mass density p 
according to 

V-[/i(||VO||/ao)VO] =4ttG P . (1) 

In the case of high acceleration, the function p = 1 and the equa- 
tion becomes the usual Poisson equation. The case of low ac- 
celeration is explained and described in a long series of papers 
(see for instance Milgrom 2002 ). At low accelerations (called the 
deep MOND regime), when ||VO|| <c ciq, the systems differ most 
from Newtonian ones: assuming //(f) ~ t, the equation becomes: 

V- [||VO||VO] = AnGa Q p. (2) 

Spherical density-potential pairs are easily built and some flat- 
tened spheroids can be numerically described (see for instance 
Ciotti et al. l20061 >. In the Appendix, we develop a new numerical 
potential solver for this equation. 

Within MOND, one finds that any disk or spheroidal den- 
sity distribution of finite mass produces a spherical potential at 
large radius (see Appendix). It is the Mondian explanation of the 
dark mater halo detected by a Newtonist. However, a disky dis- 
tribution also produces a supplementary disk potential (Milgrom 
[19831 1200T1 see also Appendix). 

Once the MOND gravitational potential is known, one can 
always apply the Newtonian Poisson equation to it, to recover 
the density distribution p„ f „ to „ corresponding to the given po- 
tential in Newtonian dynamics. Thus, how can we differenti- 
ate] Newtonian from Mondian dynamics, since the difference 
of density p newt on - Pmond (given by the two theories from the 
same potential) can always be attributed by a Newtonist to a 
dark matter component (the so-called "phantom" dark matter 
of MOND)? The answer lies within the properties of the sup- 
plementary disk of "dark matter" that appears in the phantom 

density pnewton ~ Pmond- 

The properties of this phantom disk are very specific and 
unique to MOND. Of course, two different kinds of dark mat- 
ter components (with distinct kinematics) could explain both 

2 MOND is compatible (i.e., not falsified) with recent determination 
of the Galactic halo flattening, for instance deduced from the kinematics 
of the Sagittarius stream (Ibata et al l200U Read & Moore l2005l 



the round and the thin dark components (see e.g., Read et al. 
2008 ). However, if future observations reveal such a thin dark- 
matter disk, of the exact charateristics predicted by MOND, they 
will clearly make MOND dynamics more relevant to explain- 
ing Galactic dynamics, while failing to detect it could falsify the 
Bekenstein & Milgrom ( 1984) formulation of MOND. Our goal 
hereafter is to quantify the predicted effects of this dark mat- 
ter disk on Galactic kinematics in the Milky Way. This specific 
dark (phantom) disk coud be easily identified using GAIA data; 
however, the exact predictions of MOND depend on its ability to 
describe the visible matter distribution. For that purpose, models 
such as the Besancon model are necessary because it takes into 
account important details such as the change in both radial scale 
length, and the mass distribution between thin, thick disks and 
other baryonic components. 

We note that, hereafter, when we speak about the "dynamical 
mass distribution in MOND", we always speak about the mass 
distribution that would have yielded the MOND gravitational po- 
tential from the Newtonian Poisson equation, i.e., the baryonic 
mass plus the phantom dark mass. 



3. The local dynamical surface density in MOND 

Milgrom ( 1983 ) already established that the phantom dark disk 
would enhance the K z force perpendicular to the Galactic plane. 
To be compatible with the analysis of Hipparcos data in the 
solar neighbourhood, this enhancement must not be too strong 
and must not imply an extremely massive unseen disk. Based 
on a Galactic model, Nipoti et al. (2007) determined that at 
the solar position MOND almost doubles the K z force pro- 
duced in Newtonian gravity by the baryonic disk. For this model 
(based on model 1 of Dehnen and Binney, 1998 ), the baryonic- 
disk surface-mass density is So =43 M pc~ 2 . When the disk is 
embedded in a spherical dark-matter halo, the surface density 
becomes 65 M Q pc~ 2 (see their Figure 5), while for MOND 
(baryons+phantom) it becomes 85 M pc~ 2 . This corresponds 
respectively to an increase of 51 percent and 98 percent com- 
pared to the surface density of the baryonic disk. 

Using the Milky Way MOND model of Wu et al. ( 120081 we 
found instead an increase of 57, 62, and 66 percent at R — 7.5, 8, 
and 8.5 kpc respectively (78 M Q /pc 2 at R = 8.5 kpc). This dif- 
ference with Nipoti el al. ( 12007b results from different baryonic- 
mass concentrations and different local p values between the 
models. This dynamical surface density of 78 M /pc 2 must be 
compared with the K z force determination at 1 . 1 kpc, 74 ± 6M 
pc~ 2 by Holmberg & Flynn (2004 ). Considering the various un- 
certainties, these measurements are compatible both with the 
Newton+spherical halo and MOND Besangon models. Indeed, 
these accurate determinations of the K z are based on clearly de- 
fined, but small, samples of a few hundreds of stars. Due to 
the small size of these samples, only one parameter can be ef- 
ficiently recovered: the surface density below some height z, the 
mean distance of the samples from the Galactic plane. This does 
not break the degeneracy between two effects: a small flatenning 
of the halo and a small increase of the baryonic disc density. 
These two effects can compensate without modifying the surface 
density measured at \z\ =1-1 kpc. In contrast, the K z (z) shape at 
smaller z (and also the Oort limit) will change. In the future, 
large samples from RAVE and very large ones from GAIA, will 
make it possible to recover precisely K z versus z, and to deter- 
mine dynamically the detailed vertical mass density distribution 

Pdyn(z)- 
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4. The dynamical disk scale-length in MOND 

A second and different test concerns the prediction of the 
K z force perpendicular to the Galactic plane at various 
Galactocentric radii, to measure dynamically the scale-length 
of the (visible+phantom) disk mass distribution and to compare 
this to the visible disk scale-length. The interesting property of 
this test is that only the gradients in the density distributions of 
both the visible and visible+phantom disks have to be measured. 
These gradients are certainly easier to measure than the true total 
mass distribution. 

In Fig. 1, we plot (dashed line) the radial density distribution 
of the baryonic matter in the Galactic plane for the Besancon 
model. By adopting the MOND gravitational potential (Wu et 
al. 1200 8). we then determine the baryonic +phantom dynamical 
density (solid line) that would be inferred by a Newtonist mea- 
suring the K z force perpendicular to the Galactic plane at various 
Galactocentric radii. A factor of ~ 1.25 is predicted between the 
two scale lengths, determined between 5 kpc and 10 kpc from 
the Galactic centre. 

Since the Galactic thin disc dominates the disk mass dis- 
tribution with a radial scale length of 2.5 kpc (Robin et al. 
2003), the expected disk dynamical scale-length differs signif- 
icantly, 3.1 kpc, between measurements by a Newtonist and by 
a MONDian. This measurement will be a test of prime interest 
when data from large surveys such as GAIA or JASMINE be- 
come available. 




Fig. 1. Dashed line: The radial density distribution p(R,z = 0) 
of the baryonic matter within the Besancon model. Solid line: 
The baryonic+phantom density in the Galactic plane from the 
Besancon MOND model. 



5. The tilt angle of the velocity ellipsoid in MOND 

A third new and interesting test of MOND relies on the mea- 
sure of the vertical tilt of the stellar velocity ellipsoid above the 
Galactic plane (e.g., Siebert et al. 2008). At any position above 
the Galactic plane, solar neighbourhood stars are distributed 
along isodensity ellipsoids in the (t/VW)-velocity-space (where 
U is the velocity towards the Galactic centre, V the velocity in 
the direction of Galactic rotation, and W the vertical veloicty, all 
with respect to the Sun). At a given position z above the Galactic 



plane, the inclination angle 6 of the ellipsoid minor axis with re- 
spect to the W axis (see Fig. 6 of Siebert et al. l2008b is intimately 
linked to the gravitational potential. In the f/W-plane, this angle 
is defined by 



tan2c) = 



2cr 



uw 



(3) 



where cr 2 uw , cr 2 v , and cr^ are the second-order velocity distri- 
bution moments. Cuddeford & Amendt ( 119911 1 analysed, to the 
fourth order, the consecutive moments of the Boltzmann equa- 
tion by developing these moments by expanding a Taylor series. 
Combining these moment equations, they derived expressions 
for the velocity moments and obtained an approximate expres- 
sion for the tilt angle 5, which depend only on the shape of the 
gravitational potential (but is exact only for completely separable 
axisymmetric potentials). More precisely, the ellipsoid inclina- 
tion close to the Galactic plane is related to the disk scale-length 
and also to the dark-matter halo flattening (Bienayme 2009). 
This dependence of the vertical tilt of the velocity ellipsoid on 
the shape of the potential has also been checked by numerical 
simulations (e.g., Siebert et al. 12008b . 

Here, we determine the vertical tilt angle as a function of 
Z (at the Galactocentric radius of the Sun) for the Newtonian 
Besancon model embedded in a dark halo (Robin et al. 2003), 
and compare this inclination with that obtained for the MOND 
model (Wu et al. 120081 . For this, we follow Siebert et al. (f2008 ), 
and compute the inclination of the velocity ellipsoid by integrat- 
ing orbits in the meridional plane (we note however that the 
Besancon model is non-axisymmetric, so that we average it). 
We sample the initial conditions from a Shu distribution func- 
tion (Bienayme 1 19991 1. yielding an orbit library of over 20 mil- 
lion orbits. Each orbit was integrated over 120 rotations, using a 
4th order Runge-Kutta algorithm. We then randomly selected 1 
phase-space point per orbit in the last 80 rotations. From there, 
we computed the inclination angle of the velocity ellipsoid at 
different {R,z) positions. Figure [2] displays the resulting vertical 
tilt angle (calculated with Eq. 3) as a function of z at three dif- 
ferent Galactocentric locations for both the Newtonian (+dark 
matter) and MONDian Besancon models. This tilt is also com- 
pared with the tilt expectation from a purely spherical potential 
(dotted lines on Pig . |2j : clearly, the flattening of the MOND po- 
tential due to the presence of a phantom dark disk makes the 
inclination angle lower than for a purely spherical potential, and 
lower than for the Newtonian Besangon model embedded in a 
spherical halo. The tilt differences are most significant in the in- 
ner part of the Galaxy. At the solar position, however, the dif- 
ferences remain small when \z\ < lkpc, and become significant 
only at z = 2 kpc (a difference of, 3.5 and 2 degrees, respectively 
at R-7.5 and 8.5 kpc, the inclination being smaller within the 
MOND model). At z = 1 kpc and R = 8.5 kpc (the solar posi- 
tion in the Besancon model), where the strongest observational 
constraints on the tilt have been obtained, both models predict 
an inclination of 6 degrees, which agrees with the observational 
determination of 7.3 +1.8 degrees, obtained with a RAVE sam- 
ple of 580 stars (Siebert et al. 2008). Other observational results 
on the measure of the vertical tilt remain a bit controversial and 
conflicting (Fuchs et al. 120091 Smith et al. 120091 , and will re- 
quire more (and more precise) data to avoid any potential bias, 
and confirm or disprove the sphericity of the potential. 



A measurement of the necessary accuracy will possibly be 
achieved at the solar Galactic radius with large RAVE samples 
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Table 1. Values predicted from the Besancon MOND model as seen by a Newtonist compared to observations, for the local surface 
density and the tilt of the velocity ellipsoid. 





Besancon MOND 


Observations 


E G (z =1.1 kpc) 


78 M /pc 2 


74 ± 6 M Q /pc 2 (Holmberg & Flynn 2004) 


Tilt at z = 1 kpc 


6 degrees 


7.3 ± 1.8 degrees (Siebert et al. 2008) 



in the near future, and with GAIA data in a few more years. 
It will be important to account for secondary effects (see, e.g., 
Famaey et al. 2005), such as streams, resonances, and the bar 
influence, and to correctly take into account their impact on the 
vertical tilt of the ellipsoid. However, we note, that, in the imme- 
diate solar neighbourhood, only the Hyades are clearly identified 
in the (UW) distribution that looks extremely well phase mixed, 
in contrast to the (UV) distribution where many substructures 
(from resonances or non-stationnarity) have been identified. This 
is also true for the (U, W) distributions of distant samples at 
1 kpc from the Galactic plane (Soubiran et al. 2008 and Siebert et 
al. 120081 1. The effect of the bar and specific vertical resonances 
needs to be understood in greater detail since they can modify 
significantly the vertical-tilt orientation (Olle & Pfenniger 1998). 
However, this family of orbits should also be visible as substruc- 
tures in the (UV) plane. 




-2-1012 -2-1012 

z z 



Fig. 2. Velocity ellipsoid tilts in degrees versus z at R= 6.5, 7.5 
and 8.5 (top to bottom) (the dotted line would be the inclination 
within a spherical potential). Left: within the Newtonian poten- 
tial of the Besancon Galactic model (Robin et al 2003 ), Right: 
within the MONDian potential of the Besangon model (Wu et al 
120081 . 



6. Conclusions 

In addition to the predictions made by Wu et al. ((2008 > about 
the rotation curve and the escape speed, MOND also makes 
very specific predictions about Galactic kinematics related to the 



presence of a disk of phantom dark matter. These predictions 
will be tested with future large surveys: 

(i) By measuring the K z force: at the solar radius, MOND 
predicts a 60 percent enhancement in the dynamical surface den- 
sity at 1 . 1 kpc relative to the baryonic surface density, a value not 
excluded by current data. This enhancement would become more 
apparent at large Galactic radii, where the stellar-disk mass- 
density becomes negligible. 

(ii) By determining dynamically the scale-length of the disk- 
mass density-distribution. This scale length is a factor of ~1.25 
larger than the scale-length of the visible stellar disk if MOND 
applies. This test could be applied with existing RAVE data 
(Steinmetz et al. 2006 ; Zwitter et al. 2008), but the accuracy of 
available proper motions still limits the possibility of exploring 
the gravitational forces too far from the solar neighbourhood. 

(iii) By measuring the velocity ellipsoid tilt angle within the 
meridional Galactic plane. This tilt is different for the two dy- 
namics in the inner part of the Galactic disk. However, the tilt of 
about 6 degrees at z= 1 kpc at the solar radius agrees with the de- 
termination of 7.3 + 1 .8 degrees obtained by Siebert et al. (2008). 
The difference between MOND and a Newtonian model with a 
spherical halo becomes significant at z=2 kpc. 

These tests of gravity could be applied with future GAIA 
or JASMINE data that will allow us to reconstruct the 3- 
dimensional gravitational potential of the Galaxy. To assess the 
values of the current local constraints, the predictions of the 
Besancon MOND model are compared with the relevant obser- 
vations in Table 1 . 
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Appendix A: Full multi-grid MOND potential solver 

We developed a new full multigrid (FMG) solver in the same 
spirit as Brada & Milgrom d 19951 1 and Tiret & Combes ( |2007| l. It 
solves Eq. ([T i by means of relaxation performed on the density 
field sampled at different resolutions. Large-scale feature con- 
vergence is obtained from a coarse sampling of the density, while 
high-frequency features result from fine grids (see e.g., Press et 
al. I20021 >. Relaxation is acheived with a 4 colours Gauss-Seidel, 
with 5 pre/post iterations. The solver is materially accelerated by 
means of Graphical Processor Units (GPUs) using the CUDA 
API developed by NVidia for its hardware. The Poisson equation 
is solved in ~ 1 second over a 256 3 grid using a GeForce 8800 
GTX device. Typical relative residuals obtained for the MOND 
Poisson equation are 0.4%. 

We tested our potential solver with disks of finite thickness, 
namely an exponential disk (scale length 1=1, thickness h = 0.2), 
and a thickened Kuzmin disk (scale length a=l, and thickness 
h=0.2): 



Pexp (R,z) = exp(-R/l)exp(-(z/hf) 



Pkuz(R,z) = (l+(R/a) 2 + (z/h) 2 ) 



-5/2 



(A.l) 
(A.2) 



The FMG results for the effective (visible+phantom) den- 
sity are plotted in Figs. A.l and |A.2| We note that the phan- 
tom disk clearly appears in the outskirts of the Galaxy. The 
FMG results were previously checked by solving Eq. (|2jl using 
the FreeFem++ software, a partial differential equation solver 
(Hecht et al. [2008). For spherical potentials the accuracy of the 
computed potential with FreeFem++ is excellent. With very flat- 
tened systems (axis ratio 1/10), the achieved accuracy remains 
good thanks to an adaptive grid within FreeFem++. For instance, 



in the case of the exponential disk, our mean relative accuracy 
is 1.5 per cent at radii R of between 0.4 and 6 scale lengths, 
and lower than 5 per cent for radii between 0.2 and 10 scale 
lengths. Excellent agreement is found between the FMG and 
FreeFem++, with relative differences of smaller than 1% within 
10 scale lengths. 
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Fig. A.l. Effective volume mass-density measured by a 
Newtonist from the deep-MOND potential of an exponential 
disk (scale length=2.5, scale height=0.5). 
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Fig. A.2. Effective volume mass-density as measured by a 
Newtonist from the deep-MOND potential of a thickened Kumin 
disk (scale length=2.5, scale height=0.5). 



